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I. INTRODUCTION 

Random Energy Model (REM) [1-5] is connected with the many problems of modern physics. In [6-8] has been 
found, that correlators in the directed model are connected with the free energy in directed polymer. The last is 
equivalent to REM in thermodynamic limit. 

Liouville model is closely connected with the bosonic string in the d-dimensional Euclidean space [9] . It is easy to 
check, that the connection of strings with REM is even stronger. If one considers the integration of string's partition 
via area of closed surfaces, [10-13] then after integration via zero mode of a Laplacian an expression is obtained for 
the partition like to REM with finite replica numbers (solved for real temperatures in [5]): 

Dgcj^e^ / d-^VM+QR^^ f rf^^^e"^)"^ (1) 



Here (f>{w) is a field on closed 2-d surface, a,Q are parameters real for d < 1, R is a curvature, Dgcj) is a measure. 
Q,a are defined by d according to formulas of David-Distler-Kawai (see review [12]). The analytical continuation of 
parameters Q,a at d> 1 gives complex value for parameters (see section 4). 

One can understand the last expression as an average of the fi ~ — ^-th degree (replicas number) of the sum 
e'^'f'i'^i) via normal distribution of variables = 4'{wi) with a quadratic form 
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— J d'w^g(bAq^ + QR(b (2) 

The main idea of this work (following to [6-8]) is that the phase structure of the (1) can be mapped to other models 
with the simpler choice of quadratic form in the exponent of normal distribution that the one in Eq. (2). We are 
going to connect the system (1) with the chain of models (each one with the same phase structure as the previous one 
but simpler), where the last one in the chain is the Random Energy Model (REM). That's why we decided to solve 
REM at complex temperatures [13-15] and complex replicas numbers. 
REM is a model with 
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P{E) = , exp[-— ] (3) 
> ■ 

The total distribution of energies is factorized: for 1 < a ^ /S < M 

P{E^,Ep)^P{E^)P{Ep) (4) 
The main our interest is connected with partition 

z = ^e^Y>{-(3Ei} 

i 

Z^<z^'>, (5) 

for a general value of fi. One can observe the similarity of Z defined by (5) and (1), if identify w with i, (j){w) with 
Ei, a with — /3, H = and z = J d'^w^/ge"'^^'^^ resembles J2i In (1) there is a normal distribution like (5), 

the main difference -in (1) the normal distribution is non-diagonal. In case of REM we have 2^ physical degrees 
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of freedom, like [L/aY degrees in a field theoretical model with ultraviolet a and infrared L cutoffs. The ensemble 
average (integration with a normal distribution of energies) of partition function's /i-th degree corresponds to our 
expression (1). 

In the section 2 we are going to introduce directed polymer(DP) model on the hierarchic trees with branching 
number q. The endpoints of the hierarchic tree correspond to the points Wi of the 2-d space in (1). The case g — » 1 
resembles the model (1) (for a field theoretical aspects see [17]), and the case g ^ oo is equivalent to the REM. There 
is a strict result [5] that at the case /x— > the thermodynamic limit of the introduced models are independent of q. 
In the section 3 we give a qualitative derivation of REM solution at complex temperature and replica numbers. In 
section 4 we give the classification of the phase structure of the model (1). In the appendix B. we prove, that in the 
opposite case q—^oo DP is equivalent to REM the thermodynamic limit. In Appendix A there is a rigorous solution 
of REM at complex temperatures and replica numbers. 



Let us consider the model on the hierarchic tree [2-3], [5]. Originally one has a point (origin of the tree). At the first 
level of hierarchy there are q branches. At the i-th level of hierarchy there are q new branches from the every branch 
of the i — 1-th level. At the last K-th level we have q^ end points. Let us consider field (f){x) at the endpoints. Every 
point X is connected with the origin of the tree with a single path. For the any pair of points x and x' at the level K 
it is possible introduce a hierarchic distance 



where their paths to the origin meet at the i-th level of hierarchy, V is a parameter (the maximal hierarchic distance 
between points on the tree). 

We define random variables fu on the branches at the i-th level of the tree with distribution 



II. HIERARCHIC TREES WITH CONSTANT BRANCHING NUMBER Q 



v{x,x') = 



{K-i)V 
K 



(6) 




(7) 



We define fields ^(a;) as a sum of fu along the path il{x) connecting the point x with the origin: 




(8) 



il{x) 



One can check, that 



< 'Pix^ix') 



>=V - v{x,x') 



(9) 



If one defines the distance between two points X as 




exp{v{x, x')), 



(10) 



then Eq.(lO) coincides with the ordinary expression of the 2d free field with the action (2) 



< (j){x)(l){x') >= In — 



(11) 



with ultraviolet cutoff L = cxp{V/2) and infrared one 1. 

Actually we are interested only in the distance r{x,x') for x ^ x' (for the x' = x we can take r{x,x) = 0). 
We can construct a model connected with the one defined by Eq. (1). Let us consider a partition 




(12) 



X 



At the limit ji this model has been considered rigorously in [5]. At the thermodynamic limit the model is 
equivalent to REM (with the same total number of configurations q^ and variance < E'^ >=< (/)(a;)^ >) for all values 
of q. 
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We suggest the first hypothesis of this work (it can be checked numerically), the thermodynamic limit of 
Eq.(12) is independent of parameter q lilce the case fj. ^ 0. 

We can find the pliase structure of the model at the limit q oo, it is again equivalent to REM. While the system 
(6)-(7),(12) is similar to the system (1) at every value of q due to property (9)-(ll), there is a serious difference 
between different choice of q. In the case of finite q one should consider a combinatorial problem. Only the limiting 
case g ^ 1 is similar to 2d Euclidean space, as one can use a small parameter q — 1 and construct measure. Thus we 
formulate the second hypothesis of the work: at the limit q ^ I model on hierarchic tree is equivalent to the 
system (1). 

Let us define expression (12) for the general case of q, then take the limit q ^ 1. First wc introduce the S function 
in integral representation for the partition z = exp{a^^ 4'{x)} ^ind Eq. (12) transforms into: 

/■oo 

< z'^ >= I dudv6{Rez - u)5{Imz - v){u + ivY = 

J oo 

- — 7^ / dkidk2 / dudv{u + iv)^ exp{—ikiu — ik2v)G{ki, k2) 
(27r)^ J_oo Joo 

G{ki,k2) = 

exp[iA;ii?eexp(a ^(x)) + ik2lm ex.p{a^'^ (p{x))] (13) 

X X 

Now the problem is to calculate the generating function G{ki,k2). It can be done by means of recurrence equations: 



K K 



-OO 

poo 



Let us consider the limit 



U{y) = ikiRee'xp{ay) + ik2lmexp{ay) 

G(fci,fc2) = [/if(0)]« (14) 

q ^ 1 K ^ oo 

q^ = exp{V) (15) 

Then one can express G(fci, ^2) by means of a function W{t,x): 

dW ,,,, Id^W 

dt 2 dx^ 

<t <V, -00 <x <oo 

G{ki,k2) = W{V,0) (16) 

The case of q ^ 1 trees (13)-(16) with real potential has been considered recently in [17]. Using Eq. (16) it is possible 
to found the phase structure of corresponding model in 2d Euclidean space (see. [17]). The suggested method gives 
an exact phase structure (mean field approach gives a correct list of phases but approximate borders between phases), 
as well as correct two point correlators and three point corellators for isosceles triangles [17]. In this work we are 
restricted only by phase structure. 

In principle one can find numerical solution of the last system (16) and compare it with the REM's analytical 
results for the free energy to check the first hypothesis of the work {q independence of thermodynamic limit for the 
model on a g-tree). In the next section we find the phase structure of REM, then in the appendix we prove, that 
system (14) at large q is equivalent to REM. 



III. QUALITATIVE DERIVATION OF 4 REM PHASES. 

Our goal is to calculate 

Z =< ^/'i+V2 y z = Y, 6-'^^'+'^^^^' (17) 
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where energies are distributed via (3). Let us consider these expressions for positive integer values of /i, where the 
average is over the distribution (3) for each Ei. There are two competing terms in expression of z*^ (after series 
expansion). 

The paramagnetic (PM) phase is originated from the cross terms in the series expansion expansion: 



,-dEi^-l3Ei^-..l3Ei^ 

lnZ = ,lnM + iV^=A.M±^ 

Nl3l=2lnM (18) 

The second one is the correlated paramagnetic (CPM) [4] (in Parisi's picture there is a correlation between different 
replicas), it is originated from the diagonal terms in the series expansion like to e~^^^': 

M 



Z=< (^e-^'='^')> 



lnZ = lnM+^ = Sk+^. (19) 

Let us consider continuation of (18) to the region ^ < 1. At critical temperature /3c it's entropy \nZ — 
disappears. We assume that in this region In Z is proportional to P (it is natural for a system with zero entropy) and 
fi. The continuity of In Z gives for spin-glass (SG) phase 

l7iZ = Nnl3cP (20) 

If one goes to complex temperatures [11-12], then (18) transforms to (it is easy check directly for integer n) 

,^Z = N^-^±^Mt (21) 

For the SG phase one has to replace /3 by (3\ in (20): 

InZ = N^J^|3c|3l (22) 

For complex temperatures there is a fourth, Lee- Yang-Fisher (LYF) phase. The derivation is not direct. The point 
is, that for noninteger values of ^, 

Z'^<\zY-> (23) 
After this trick it is easy to derive the LYF expression. The principal terms are e~^'^i^': 

lnZ=iM±M)^ (24) 
4 

Let us now continue our four expressions to complex values of [i. For PM phase an analytical continuation of Eq. 
(18) gives 

l^^^^(/3g+/3?-/3|)Mi-2/?r/3.M2 (25) 



For SG phase we have 
For LYF phase: 



lnZ = Ar/xi/3c/3i (26) 



lnZ=^M±M)^ (27) 
4 



For CPM an analytical continuation of Eq. (19) gives 



^^Z= + {0i - P2M - Mi) - 4/3i/?2MiM2)] (28) 
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To find the borders between four phases one should first find the correct phase at /i ^ limit, then compare its finite 
/Lt expression for | In Z| with the corresponding one given by CPM phase. It is known, that LYF phase exists at [14,15]: 

(29) 

and PM one at /3 < /3c. For a complex temperatures one has a condition for SG phase 

/3i > /?c + 01. (30) 
The last point. Strict derivation gives, that LYF for noninteger ^\ exists only at 

Ml > -2 (31) 
The paramagnetic phase is the most symmetric one, there are local symmetries in the model. 

For the case of SG phase there is some order and no local symmetry. For the c;ase of LYF phase there is some 
correlation between couples of replicas. Those two phases (SG and LYF) resemble non-unitary models in field theory. 
For the case of CPM phase there is a correlation between all the replicas, but local symmetries are conserved. This 
phase has not any pathology. Thus together with PM phase it can be connected with unitary models. 



IV. STRING PHASES 



What can we say about quantum field theory and strings on the basis of our results? One should understand that, 
in some sense, there is a hierarchy of requriements physical theory to be mathematically rigorous. 

On the top level in quantum field theory one demands the unitarity of the theory as the main constraint. In the 
lower level of statistical field theory in Euclidean space the genesis of two probabilities: in ensemble and Boltzmann one 
is not too complicated. Therefore one can consider non-unitary model also, connected with some physical situations 
[18]. We have another physical constraint: finite number of primary fields. In the lowest level of hierarchy is an 
ordinary statistical mechanics. Here one is happy, when can construct a thermodynamic limit {N — > oo) for the free 
energy, entropy. 

We see, that our results are quite restrictive. First, we should forbid a situation like Lee- Yang-Fisher singularity 

with too negative replica numbers /ii < —2. Here it is impossible construct a thermodynamic limit. Other situations 
with Lee- Yang-Fisher phases as well as with spin-glass are a bit interesting, but here could not be any unitary theory. 
The most interesting are paramagnetic and correlated paramagnetic phases. In this area there is a some chance for 
unitarity. 

Let us return to the partition of bosonic d-dimensional string (1). For the ultraviolet cutoff L and infrared one a 
the number of degrees is 

7-2 

M=— . (32) 

a, 

configurations. Let us define distribution of (l){w) over all points w, using the free field action from (1): 

dp' 

p{M =< S{<f>o - <f>{w) >^(y,)^ ^""^^"25^^' ^^^^ 
where G is correlator of (f>{w) fields, the average is over the distribution 

p((/)(w)) ^e^-f d'^Vg'PA>t>+QR<i>^ (34^ 

and 

G(0)=21n^. (35) 

We replace our system (1) with a REM model having the same number M independent variables Ei ~ 4>{w) with the 
same distribution (32): 



lnZpM = 2^4^1n- (36) 
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Here Zpm is the partition connected with the Eq. (1) for the PM phase at real a. We rescale the temperature: 



^=/3,/3e = l,M = -- (37) 



DDK formulas give for d = c dimensions: 



Q=\l^-,a = --={V25^c-VT^c) (38) 



For the sphere topology: For the 1 < cJ < 25: 



,P2 = - 



24 V24 



Ml = Y^[25 - c],M2 = V(25-c)(c-l)- (39) 



For 25 < < 26 we have: 



^V^£25_x/^ 
\/24 

Ml = ]^ [25 - c] + V(25-c)(c-l)^ (40) 

For other string topologies one should rescale the /x expressions in Eqs. (39)- (40): 

^,^{l-g)^l (41) 

Let us consider first the sphere topology. We denote y = ^fj^- For the 1 < d < 25 we have: 
Pi = ^Jy,(i2 = — \/l — y, Ml = 2y, M2 = 2-\/y(l — y). We derfive the following 4 expressions for the InZ: 

Mi[l + (/32 - PI)] - 2H2P1P2 4y2 + 4y(l - y) 



2 2 
1 + (M? - Mi)(/3f - - 4miM2/3i/?2 _ 1 + 42/(2t/ - 1)^ + 16y^{l - y) _ 1 + 4y 



2y,PM 
,CPM 



2 2 2 

Ml/3l = 2y3/^5G 

iMl±^^yil±M,LYF (42) 
4 2 ' ^ ^ 

We see, that CPM phase is prefferable in the region 1 < < 19 (0 < y < 1). 

For the 19 < d < 26 we should compare \nZ expressions for the CPM and LYF phases, as /3i = 0,mi > 0. Now 
we denotee y = -^|^,mi = -2y + 2i/y(y + 1),/32 = y - \/rTy. For the CPM phase we have InZ = = 
(1 - 4(y - ^y(y + l))2(j/ - and for the LYF phase \nZ = if = -'^v+^VyT^ _ ^e see that for spherical 

topology for the whole region 1 < d < 26 string is in CPM phase. 

Let us consider the torus topolgy case. We have for PM phase InZ = — ^'"^^ = y, for SG phase In Z — Pi = ^/y 

and for LYF phase \nZ = = where y ^ At 1 < d < 19 system with torus topology is in the SG 

phase. LYF phase exists at 19 < d < 26. 

Let us consider now now higher topolgies g <2. Again system is in SG phase for 1 < d < 19. At 19 < d < 25 still 
exists a thermodynamic limit and system is in LYF phase, if ^^f^{g — 1) < 2, therefore g = 5, d = 19 is a multicritic 
point. For the d = 26 there is a thermodynamic limit with LYF phase at {g — l)/3 < 2. 

Let us consider now a case of superstring. Now one has: 



Let us denote u = 



19 -d V9- d - zVc^ 
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Ml = 



9-d 

4 



= 2m, /i2 = 



V(9-rf)(l-rf) 



2^/u{^ - u) 
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We see a mapping y ^ u. Now the transition to the LYF phase is at d = 7. According the [19] interesting dimension 
is = 5, connected with QCD interpretation as strings [20]. 

What one can say about string's physics on the ground of the REM picture? The most interesting case is the sphere 
case. When one chmbs over the d = 1 barrier, nothing happens in REM picture, system is still in CPM phase, as for 
the d < 1. The free energy has not any singularity (might be there are singularities in some correlators). To reveal 
interesting (unitary) theories explicitly one should solve the directed polymer at finite replica number including finite 
size corrections and correlators. But at least for the sphere case the REM analysis seems to be quite reliable. 



In sections 1,2 and Appendix B we gave an arguments for the connection of string's partition with a finite replica 
number REM. In section 3 and in Appendix A we solved Random Energy Modela at complex temperatures and replica 
numbers. In section 4 we take string model with an analytical continuation of the David-Kawai-Distlcr formulas at 
d> 1 and mapped it to REM. The validity of DDK formulas at d > 1 is still under question, but we hope that the 
analytical continuation could reveal singularities of the system. It is a typical situation in statistical physics, when 
there is a singularity in free energy expression, when it is analj^tically continued from one of phases to the border 
between phases. 

We obtained a bit strange result about difference of phases for the different topologies of string surfaces at d > 1. 
For the spherical case there is no any barrier at d = 1, at least for the free energy. For the other topolgies at d > 1 
system is in SG or LYF phase, sometimes the model is so pathalogical thet there is no any thermodynamic limit. 

There have been early attempts to connect strings with spin glasses. I have several discussions with V. Knizhnik in 
Alma-Ata conference in 1985, later in Yerevan before his death. He was highly intrigued with a ultrametry property 
of spin-glasses and trying to connect them with string. 

M. Virassoro also informed me about his and G. Parisi's attempts to connect strings with spin glasses. 

We could succeed due to work done in [4] , [6-8] and a simple observation that string's partition is similar to finite 
replica REM just after zero mode integration. String theory is too mathematized. In this work we tried to catch some 
narrow but crucial aspect of the theory using more physics and less complicated mathematical tools. 

I am grateful to B.Derrida, Y. Sinai for the discussion of g = 1 trees. 
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V. 



CONCLUSIONS 



APPENDIX A: REM'S SOLUTION FOR COMPLEX TEMPERATURES AND REAL REPLICA 

NUMBER. 

To calculate expression (5) we introduce an identity 

J dUiSiUi - ReY,e^^'+'^^^^' j dU^SiU^ - ImY^e''^^+'^'^^* = 1 

i i 

and an integral representation for S function 5{z — u) = J rffce*'^^'^""^: 

f{h,k2)=9{h,k2)" = 

^= / dxexp[—]exp{ikief^^''cos{P2x)+ik2e^'''sin{p2x))]'^ (A.l) 
V Nit J -oo ^N 

Having an expression for the function /(fcl, k2) we can define the partition Z: 

Z = 

dkidk2dUidU2e-^^'^^-^^^^^{Ui+iU2Yf{k^,k2) (A.2) 

J — oo 

This is an exact expression. In thermodynamic limit we will consider four different asymptotics for the function 
/(fcl,fc2). 

In the paramagnetic phase we expand an exponent via degrees of fci, 

g(fci , fe) « 1 + ikiRee^ ^ + ik2lme^ 5 (A.3) 

Integration via dki , dk2 gives 

5{Ui - i?ee^^^^ )S{U2 - Ime''- 

Eventually we derive for the PM phase: 



f{ki, k2) « exp[ifciMi?ee^ s + ik2MIme^ 

We miss the imaginary part in the expression of In < Z > . 

For the Lee- Yang-Fisher (LYF) phase we take the second terms in the expansion of the exponent : 

^(fc ) ^ 1 3_ r ^^+2/3.x (fcl C0S{(32X) + k2Sin{(i2X)f 

= 1 - ^l±^e2^^i (A.5) 



We obtain: 



/(fci,fc2)«exp[-M^^l±Me2^/5?l 



Z 

2 I Z-2 



^ i.2 I j,2 

/ dA:idfc2dt/idC/2e-*'=i^i-^'=^^^(C/i + ^C/2)^'^l±^ exphMe^^'^' 

= '^f^^^^^exp[-M±^](t/,+^f/2r-- 
= e^iN^Hr^_ dy)exp[-r2]r'^i+i 

TT Jo Jo 



_ l.^^,jv£l+lg£p. Mi + l + ^M2 . exp(2^(/ii +i/x2))- 1 

TT ^ 2 ' IJ.i+ilJ,2 
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The LYF phase (A5) is rigorously defined for /ii > —2, otherwise there is a singularity at integration dUidU2- 

For the CPM phase we consider the case of positive /Xi. We take n>/Xi>n— l,i^ = /Ui — n, 0>i^>— 1 and write 
an equivalent expression for the (Al): 

J-oo J-oo dki dk2 

{Ui + iU2r (A.7) 

Now we assume, that in the principal region of dki, dk2 integration 

/(fci,fc2)-l< 1. (A.8) 
Therefore we can expand the exponent in the f{ki, k2) expression: 



1 [°° 

/(fci, ^2) = [ / da;exp[-— + ifci cos(/?2a;) + ik2sin{l32x)\[ 



M 



-^J-dxe^p[-^ + ie^'^'iki cos(/32a;) + k2Sinif32x))] - 1} 
-00 VttN Af 

Then after integration by parts: 



M r 



dk^dk2dUidU2e-'^^'''-'^''-{Ui+iU2Y{-^^ + -^TfikiM) 

M f°° f°° 1 

= — r/ dkidk2 dUidU2e-'^^^^-'^^HUi + iUiY ^=dx 

exp[-— + iRe{ki - iA:2)e(^i+*^^)"]e('^i+''^=)^" (A.9) 

We miss the term 1 in the expression of f{ki, k2), because its contribution is equal to after integration by parts. Let 
us denote E = exp[(/3i + il32)x],K = kexp{iip) — ki + ik2, U = Ui + iU2- First we take the integration via dki,dk2- 
The result is 6{E — U). Then we calculate Gaussian integral via dx and derive an expression for the correlated 
paramagnetic phase (CPM): 

1 r°° 2 

Z =< z'' >= M^= / rfa;e-7r+/'^(/3i+»/52)^ 



,,pt^(MrMi)(/3r/3|)-W2MrM.+/3gj (A.10) 

Let us calculate SG phase. It is convinent to use another representation of function f{ki,k2) [16]. Using the 
Stratanovich transformation for the energy density term 

exp(-— ) = -j=r- / dye ^ +/*iV^xy 

^ V^TT J-ioo 

we derive 

foo /.oo ATo2.,2 



/AT /-OO /•OO AT OZ Z 

g{ki,k2) = ^ / dy dxeM—^ + PiVNxy + ■ikiRee^^'+'^'^^'' + ik2lme^^'+'^'^'^^) (A.U) 
After transformation v = k exp{-</N pix) we have 

fih,k2) = {— dye-y'^''+^G{y,k,^)) 

POO 

G{y,k,ip)= (^^e™cos[/32//3i(ln^-lnfe)-v.]^(y-l)_ (^_;^2) 

Jo 
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Wc are interesting in Eq.(12) for the | lnfc| ^ A'', therefore we can calculate the asymptotic of the function /(fcl, fc2) 
via the saddle point method. There is a pole of function G{y) at y = with residue equal to 1 (it can be derived 
putting a smal low integration limit). Let us shift the integration loop to the saddle point. We have: 

-I pioo ^2 2 

/(fc,^) = M(l + — / dj/e-^i"l'=le^^G(2/,A:,^)) (A.13) 

For the saddle point we have 

.0 = ^ (A.14) 

We move the integration loop via dy to catch the saddle point. For the analytical continuation to the region [— 1 < 
Rey < 0] we transform expression of G{y) from (A12) using the integration by parts: 

G{y, k, <p) = -i dvvy exp[it; cos [/^s//?! (In - - <^]]{cos[/J2//3i(ln ^ - <p] - sin[/32//3i(ln --^]} (A.15) 

We have an asymptotics: 

1 9 In t inii^ 

/(A:i, k2) = exp[-Me A] 

A = -G{yo,k,ip) (A.16) 

One should take only this asymptotic instead of (A3),(A5) if, while shifting the integration loop, we don't intersect 
the pole at y = — 1, 

"^"^"<1 (A.17) 



Nf3f 

Otherwise, at ^^0- > 1 we should consider all three different asymptotics (A3),(A5),(A16) and choose the largest 
one. From the Eq. (al6) we derive immediatly the bulk value of partition: 

Z ~ exp[MiAr/3i/3e] (A.18) 

We derived accurate expressions for the PM phase (A4),LYF phase (A6), CPM phase (AlO) and bulk expression for 
the SG phase (A17). We see, that LYF phase can be constructed only at iii > —2. Thus a situation, when bulk 
expression for the < Z > is given by LYF phase and /zi < —2, model is to pathalogic and there is no thermodynamic 
limit. 

To find the borders between phases one should solve the model at the limit jj, ^ 0, then compare the corresponding 
expression of the largest free energy with the one in CPM for the finite (it exists only at fj,i > 0). One should 
choose a phase, having larger value of | lnZ|. Let us remember also that SG phase exists at /3i + /?2 > Pc, the LYF 
phase at /32 > , /?! < ^ • 

APPENDIX B: EQUIVALENCE OF REM AND DIRECTED POLYMER IN 

THERMODYNAMIC LIMIT 

We consider a hierarchic tree with K levels and large Q. We have, that along any path connecting endpoint with 
the origin: 

J2<el>=N (B.l) 

a 

Let us first consider the PM phase. We should define the generating function 

< f{ki,k2) >= 

< exp[iiie(fci - ifc2)e''i+'''=)(^'i+^'2+-«i^.)] > (B.2) 
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Here the sum is over all the paths, connecting cndpoints with the origin. Let us consider first the integration via the 
last level of hierarchy. We expand an exponent and after integration via dej^ : 

< G{kxM) >=< exp[ii?e(A;i - ik2)e^'^'+''^'^^''+'^+-"'-'^]Qe^^^^^i^^^ > (B.3) 
Repeating this procedure K times, we obtain 

< /(fci, fe) >= exp[iRe{ki - ik2)Q^e^ ^ ] (B.4) 



In principle the expression in the exponent could be large. We recover the REM result for the PM phase with accuracy 

0(1). 

For the LYF the integration via the last level of hierarchy gives 

< /(fci, fc2) >=< 1 - ^L±Me2ft(ei+^2+...^-i)QM±Me^] > (B.5) 

Repeating the integration K times (we expand the exponent all the time via {kf + fc^)) gives the REM result for LYF 
phase 

< /(fci,fe) >= exph^^^Q^e^/'i] (B.6) 
Again we have equivalence with accuracy 0(1). 

The case of SG phase is a bit complicated. Now the integration via the last level of hierarchy gives 



< f{ki,k2) >= exp[- 



^ [lnfc + /3i(ei+e2 + ..tg_l)l'' 



Nc 



-e 



k = Jkf + kl (B.7) 



where c 0{\). This expression resembles the case of real temperatures, where the generating function is calculated 
[16]. Using those results for the real temperature directed polymer, we derive again the REM expression for the SG 
phase. 

Let us consider now the case of GPM phase. Here we use again formulas (a7), (a8). Now all the integrations 
decouple and we recover the result of a simple REM: 

< 2^ >= M < e~(''i+'''2)('^i+*''2)^ >= 
cxp[.V ^^' ~ ~ ~ 4/?i/?2MiM2 



